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We review efforts to study, using the methods of lattice field theory, the phe- 
nomenon of diquark condensation via BCS pairing at a Fermi surface, which has 
been proposed as a mechanism for color superconductivity in dense quark matter. 
The particular models studied are the Gross-Neveu model and SU(2) lattice gauge 
theory; in both cases evidence for superfluidity at high density is presented. The 
behaviour expected for quarks in both fundamental and adjoint representations of 
SU(2) is contrasted. 

1 Introduction 

At low temperatures, fermionic matter becomes degenerate; ie. all single par- 
ticle states are occupied up to some sharply defined scale, the Fermi energy, 
which characterises allowed physical processes. In a normal Fermi liquid, ex- 
citations above the Fermi energy can have arbitrarily small energy. However, 
if there is an attractive interaction between fermion pairs, and the Fermi en- 
ergy is large enough, then the normal state should be unstable with respect 
to the opening of a gap at the Fermi surface. In field-theoretic language, this 
gap is signalled by a non-vanishing condensate formed from two elementary 
fermion fields. If our fermions are quarks, then we have a diquark condensate 
(qq) 7^ 0. As reviewed below, diquark condensation can lead to new and more 
interesting ground states. In this talk we will summarise our efforts in, and fu- 
ture prospects for, probing this phenomenon using the methods of lattice field 
theory. Throughout we will be attempting to explore the zero temperature 
limit, although this may be difficult to reach in practice on finite volumes. 

When considering the possible diquark condensates that may form, vari- 
ous questions present themselves. Firstly, if the quarks carry a gauge charge, 
is the condensate gauge invariant? If so, then the state with (qq) ^ breaks 
some global symmetry of the Lagrangian, implying eg. that excitations carry 
indefinite baryon number. This is the phenomenon of superfluidity, which 
occurs in condensed matter physics at milliKelvin temperatures for He 3 . The 
systems we have studied to date are expected to exhibit superfluidity. How- 
ever, if the condensate breaks a local symmetry, then the system will be a 
superconductor, and the diquark condensate a direct analogue of the Cooper 
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pairing found in metals at low temperature. Since the quarks of QCD trans- 
form under a non-Abelian symmetry, the phenomenon is still more exotic: 
color superconductivity has been postulated as a dynamical embodiment of 
the Higgs_mechanism, making some or all of the gluons massive at high baryon 
densityJiltrcl Secondly, is the condensate a spacetime scalar? Naively one might 
expect so, but examples of rotationally non-invariant condensates are knawn 
in condensed matter physics (He 3 A), and have been postulated in QCD.I3 A 
parity violating diquark condensate is another possibility. 

The most crucial consideration is that the condensate must respect the 
Pauli Exclusion Principle, so that its wavefunction is antisymmetric under 
exchange of all available quantum numbers, which in the case of quarks are 
spacetime, color and flavor. This has the consequence that the ground state 
may be extremely sensitive to the number of light flavors present ,□ or as 
we shall see below, to the representation of the color group carried by the 
quarks. This impact of the Exclusion Principle makes the phenomenon of 
intrinsic theoretical interest, and may even underlie the known difficulties in 
simulating QCD with non-zero baryon density.cl 

2 Four Fermi Models 

The first model we will consider is referred to as either the "Gross-Neveu" 
(GN) model in 2+1, or the "Nambu - Jona-Lasinio" (NJL) model in 3+1 
dimensions, and is a relativistic generalisation of the BCS model originally 
used to describe superconductivity. It has the following Lagrangian density 
in continuum notation: 

C = V# + m)if> - g 2 [(^) 2 - (^75^) 2 ] , (1) 

where ip is an isodoublet transforming under a global axial SU(2) L ® SU(2) fl 
symmetry, which is exact in the chiral limit m — > 0: 

ipL 1 * UipL , ipL 1 y 4>lU^ ; ip R i-> Vip R , ipR h-> iJ)rV*. (2) 

C is also invariant under a global U(l)y of baryon number: 

; ^^^g-^. (3) 

Let us review some properties of the model. First, for sufficiently strong 
coupling g 2 , the axial symmetry spontaneously breaks to a S\J(2) V of isospin 
by the formation of a chiral condensate (tptp). The spectrum in the chirally 
broken phase contains massive "baryons" - the elementary fermions, and iptjj 
composite "mesons" , which include three light Goldstone pions. The chiral 
condensate sets the dynamical mass scale M/A for the model, which can thus 
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be taken to zero at some critical g 2 , denning a continuum limit. A special 
feature of the GN model is that the continuum limit defines an interacting 
field theory, whereas that of the NJL model is logarithmically trivial. a For 
the present purposes, the most interesting feature of the models is that they 
can be formulated on a lattice, and simulated with baryon chemical potential 
fx 7^ by standard Monte_Carlo methods the reasons for this are subtle and 
only recently understood.Q It is found that for /i greater than some fi c of order 
the baryon massJhe symmetry (||) is restored. 

In our worlo3 we have focussed on the formation at large \i (for non- 
interacting fermions /i is precisely the Fermi energy) of a diquark condensate 
(qq), where in continuum notation the diquark wavefunction in 2+1 dimen- 
sions reads 

qq = iJj tr (Cj 5 )<S>T2<S)T2ip. (4) 

Here C is the charge conjugation matrix, and the C75 structure ensures that 
the condensate is a scalar. The first of the T2 matrices operates on an im- 
plict flavor space due to the species doubling generic to lattice fermions, and 
the second on the explicit isospin indices. Since all three operators are anti- 
symmetric matrices, the overall antisymmetry of the wavefunction is ensured. 
This condensate respects the axial symmetry (Q) but spontaneously breaks 

u(iVi). 

2.1 Two-Point Function Approach 

Our first attempt to look for a signal for (qq) 7^ was via the asymptotic 
behaviour of the diquark propagator, which by the cluster property should be 
proportional to the square of the condensate: 

G(x) = (qq(0)qq(x)) = (qq(0)qq(x)) c + (qq)(qq) => km G(x) = \(qq)\ 2 . (5) 

x — >oo 

The condensate should thus manifest itself .as a plateau in the large-f be- 
haviour of the timeslice propagator. Resultsa for G(t) — ^2$G(x,t) from a 
16 2 x 40 lattice, with [i = 0.8 > n C) are shown in Fig. 1, clearly showing a sta- 
ble plateau when compared to the same quantity obtained for non-interacting 
fermions, shown with closed symbols. Taking the square root of the plateau 
height as a measure for (qq) results in the plot shown in Fig. 2, where chiral 
symmetry restoration at /i c ~ 0.65 followed by a rapid rise in baryon number 
density is also clearly visible. We see a dramatic increase in (qq) proceed- 
ing from the low density chirally broken phase to the high density chirally 
symmetric phase, suggesting that diquark condensation is taking place across 
the transition, as revealed by long range order in the timelike direction. This 
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Figure 2. 



result is consistent with the notion of competition between chiral and diquark 
condensates] 2 ! We have also observed signals both for a small scalar diquark 
condensate in the broken phase, and for a non-vanishing pseudoscalar con- 
densate in the dense phase, implying spontaneous breaking of parity; both 
these may well be finite volume artifacts, however. 

Unfortunately, data from lattices with spatial volume L? s varying from 
8 2 to 24 2 do not support a naive application of Eq. (||) , which suggests the 
plateau height should be an extensive quantity. Indeed, the plateau height 
saturates for L s > 16, a result which may be due to the influence of Goldstone 
fluctuations, which would wash out the signal in the absence of an external 
source, or from the relatively small number of participating qq states in the 
vicinity of the Fermi surface on these moderate systems. 

2.2 One-Point Function Approach 

In an attempt to clarify matters we have also performed direct measurements 
of the diquark one-point function (qq). By rewriting the fermionic action 
Sferm = t^Mip in the Gor'kov representation, it is possible to add explicit 
diquark source terms: 

Z[j,J] = (Pt4[j,j|>. (6) 
The diquark condensate is then defined by 

^-AHH^i 0})- (7) 



(qq) 
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Figure 3. Figure 4. 

which is straightforward to implement. Our results for (qq(j)) in the GN 
model from a 16 2 x 24 lattice are shown in Fig. 3. We see a large jump 
in the signal as the chiral transition takes place; at low density the signal is 
approximately linear in j, whereas at high density there is a marked curvature. 
In neither case, however, is there convincing evidence for a non-zero intercept 
as j — > 0. 

It should be stressed that these measurements are "quenched" in the 
sense that the diquark source term is not included in the Monte Carlo update 
algorithm. We are currently attempting full simulations with j J= in the 
expectation that finite volume effects will be easier to interpret .0 Inclusion 
of the source term will also enable estimates of disconnected contributions to 
the two-point function, which may have a better overlap with the Goldstone 
mode if one exists. It will also be valuable to develop spectroscopy techniques 
in the Gor'kov representation, since in the presence of a gap the eigenstates of 
the transfer matrix will be a superposition of q and q states. It is safe to say 
that much more work is needed before the apparently straightforward exercise 
presented by four fermi models is complete. 

3 Two Colors QCD 

3.1 Fundamental Quarks 

SU(2) gauge theory with staggered lattice fermions in the fundamental repre- 
sentation does not suffer from the difficulties associated with simulatingjdense 
QCD, since it can be shown that detM is real^ and positive for all /Jill In 
Fig. 4 we show results from runs on a 4 4 lattice, with Nj — 4 physical flavors 

a This follows readily from detM = detT2Mr 2 = detM*. 
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of fermion and m = 0.2, in the strong gauge coupling limit as a function of \x. 
For n < 0.4, the system remains essentially unchanged, but for 0.4 < /i < 0.9 
the chiral condensate (^i/j) smoothly decreases to zero, at the same time as 
the baryon density n = (tpjot/j) increases from zero to its saturation value of 
two quarks per lattice site. The numerical effort required, as measured by 
the number of computer iterations required to invert M, rises steeply in the 
crossover region. Most interestingly, the average plaquette decreases as den- 
sity rises, until at saturation it assumes the expected quenched value of zero. 
This is a signal of Pauli blocking; at high density virtual qq pairs are sup- 
pressed for kinematical reasons, and color screening via vacuum polarisation 
thus reduced. 

A physically appealing way of understanding the rolexif the complex phase 
of detM in gauge theories is in terms of conjugate quarksS3 Monte Carlo simu- 
lations of QCD demand a real functional measure, and therefore update using 
det MM t = detMdetM*. The M describes quarks q in the 3 representation 
of SU(3), the M* conjugate quarks q c in the 3. Usually conjugate quarks 
are regarded as anti-quarks of a different flavor, but once jj, ^ it can be 
checked that both q and q c couple to [i in the same way and hence carry pos- 
itive baryon number. Therefore simulations with a real measure can describe 
gauge singlet qq° bouncLstates; in particular a pseudo-Goldstone state, the 
baryonic pion can formjla and is the lightest baryon in the spectrum, since by 
the usual PCAC arguments its mass is expected to vanish as y/m in the chiral 
limit. Now, simple energetic arguments suggest that observables should start 
to show /x-dependence for /i > /i G ~ mib, where mib is the lightest baryon 
mass. The presence of an unphysical light baryon therefore plays havoc with 
the physics of ^ ^ 0.[] 

For SU(2) gauge theory, however, baryonic pions are not unphysical. Con- 
sider the kinetic term for staggered lattice fermions: 

Sun = 2 V/m(x) [x{x)U ^{x)x{x + (i) - x{x)Ul{x - fi)x{x - A)] • (8) 

0,3 

We can identify both U(l)y and U(l)yt global symmetries, which may be 
associated with baryon number and axial charge conservation respectively: 

U(l)y : X >-> e M X X ^ Xe~ la U(1)a ■ X -> e^ X X ~ X^ e , (9) 
where e{x) = (— l) x o+ x i+ x 2+ x 3. However, by defining new fields 

X e = { Xe ,~X t lT2), = ( X t o r ,XoT 2 ), (10) 

b For four fcrmi models, qq states are much lighter than qq c due to disconnected diagrams. @ 
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where the subscripts distinguish between even and odd lattice sites, then due 
to the pseudoreality of the 2 representation the action may be recast as 

Skin = x X! ^^ X ) [X e {x)U IJ ,(x)X (x + p) - X e (x)U^(x - p)X (x - p)] . 

xeven.fi 

(ii) 

Since the X and X fields have two components, the global symmetry 
is enlarged to U(2), and can be viewed as relating mesonic qq states to 
baryonic qq states. In particular, a condensate (xx), which in the chiral 
limit at zero density would be expected to break spontaneously the orig- 
inal U(l)^4 symmetry, can be U(2)-rotated into a gauge invariant diquark 
condensate (qq2) = (x trT 2X), which also breaks U(l)y. In cither case the 
overall symmetry breaking is U(2)— >U(1), which is different from the_pattern 
SU(2N f X.-> Sp(2iV/) predicted for continuum SU(2) gauge theory.tl A full 
analysis&j shows that in the limit m — > 0, — > the symmetry breaking is 
accompanied by three Goldstone modes, with quantum numbers 0~, + , + . 
Depending on the direction on the U(2) manifold chosen by the condensate, 
the Goldstones can be regarded as either mesons or baryons, and for this 
reason the presence of baryonic pions in the Monte Carlo simulation is not 
harmful.^] 

Once m and /j, differ from zero, the high degree of symmetry is no longer 
present. Let us first present some simulation results obtained using a hybrid 
Monte Carlo (HMC) algorithm, which corresponds to Nf = 8 physical flavors. 
In Fig. 5 we show results for (qq2(j)) obtained using eq. (Q) for various fiona 
6 4 lattice with quark mass m = 0.05, and gauge coupling = 1.5.E3 There is a 

c For the model described by the measure detMMT the actual pattern is U(4)— >0(4). 
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clear distinction between the curves for fj, < 0.4, which appear to extrapolate 
to zero as j — > 0, and those for fi > 0.5, which plausibly yield (qq) ^ in 
this limit. Results for (qq2) obtained using a simple polynomial extrapolation 
to the zero source limit are plotted as a function of fi in Fig. 6, along with 
the chiral condensate (xx) an d the average plaquette. The chiral condensate 
appears to decrease smoothly as soon as \i > 0, until by \i = 0.8 it actually falls 
below its free-field value. The immediate onset of the fall with /i is difficult to 
understand unless there are thermal effects associated with the rather small 
lattice volume. The diquark condensate, by contrast, stays small (probably 
even zero within systematic errors) until a value /i c ~ 0.4, whereupon it 
rises sharply to a large non-zero value, and then almost immediately starts 
to fall again. Evidence for a phase transition at /i = \i c also comes from the 
average plaquette, which shows a discernible kink between the low density 
phase, where it is approximately constant, and the high density phase where 
it decreases with /i, presumably as a consequence of Pauli blocking. 

It is tempting to associate the decrease in (gefe) with /i in the dense phase 
with asymptotic freedom, since the scale defined by the Fermi energy rises 
monotonically with /i, and further measurements reveal that n also rises sig- 
nificantly over this range of fi. It will require much more extensive simulation 
on a range of lattice volumes and spacings, however, before this effect can be 
distinguished from lattice artifacts due to saturation. 

As mentioned above, the diquark condensate spontaneously breaks the 
original U(l)y symmetry of @, which remains a symmetry of the Lagrangian 
even once m, /J, ^ 0. Therefore we expect the appearance of the condensate 
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in the dense phase to be accompanied by a scalar Goldstone mode. In Fig. 7 
we plot the contributions to the diquark susceptibilities, 

X« = X>«(0)M(a!)>, (12) 

X 

from connected quark line diagrams in both scalar and pseudoscalar channels. 
Both increase sharply for \l > /i c ; the scalar signal, which should project onto 
an exact Goldstone mode, exceeds the pseudoscalar, which projects onto a 
pseudo-Goldstone mode. The two channels should become degenerate in the 
chiral limit O Once again, it is likely that to establish the level ordering in 
the dense phase unambiguously will_require the contribution of disconnected 
diagrams to be taken into account £3 

Fig. 8 is a suggested phase diagram for SU(2) lattice gauge theory with 
fundamental quarks in the (to, \i) plane. In principle there are three phases: 
the vacuum, for which both baryon number density n and diquark condensate 
(qq) vanish; a normal phase for which n > but (qq) = 0; and a superfluid 
phase^] (recall qqi is gauge invariant) for which both n > and the order pa- 
rameter (qq) 0. Note that the chiral condensate (xx)> which is usually con- 
sidered as an order parameter, is everywhere non-zero in this plane, although 
from Fig. 6 we expect it to decrease from bottom right to top left. The dashed 
line separating the vacuum from the normal phase is /z (m) oc a/to, following 
the arguments about the lightest baryon given above. The solid line iL c (m) 
which separates the superfluid from the normal phase may well coinckLe with 
\x for some or all of its extent. It will be a goal of future simulationstS to es- 
tablish in the first instance whether the normal phase exists, and if so whether 
it is confined to the large m, [i region where (xx) is favoured kinematically, 
but (qq) suppressed by asymptotic freedom, or extends in a narrow tongue all 
the way to the zero density chiral limit at the origin. Another issue to explore 
will be the persistence of the superfluid phase for T > 0. Finally, even though 
the model has the wrong physics for fi < u c , studiea of gluonic dynamics in 
the dense phase may still be of relevance for QCDJIZI in particular we might 
expect competition between the enhanced screening expected from a non-zero 
density of color sources in the ground state, and the kinematic anti-screening 
due to Pauli blocking. 

3.2 Adjoint Quarks 

The consequences of formulating the theory with adjoint rather than funda- 
mental quarks are profound. Even at zero density, the adjoint model is distinct 

d Unlike He 3 , superfluidity arises here via BCS pairing in the s-wave channel. 



trento: submitted to World Scientific on February 1, 2008 9 



because of the possibility of gauge invariant spin-i bound states, either qg or 
qqq, in the spectrum. We wish to advocate SU(2) lattice gauge theory with 
adjoint quarks as a "Toy QCD" for the purposes of non-zero density studies. 
First let us discuss the symmetries of the kinetic term for a single stag- 



gered flavor. The manipulations leading to eq. (11) go through as before, but 



now the can be chosen real, and the X, X fields defined: 

X e ^(xe,xT), Xl r = (xl r ,Xo)- (13) 

Formation of a chiral condensate (xx) causes a U(2) — > Sp(2) symmetry 
breaking, resulting in one broken generator whose Goldstone can be identified 
with th£ familiar pion. Once again, this is the opposite of the continuum 
result Ilia Since there are no light diquark states in this case we expect no early 
onset, ie. lim m ^o Mo 7^ 0. 

Now consider possible diquark condensates which might form at high 
density. In the absence of a detailed dynamical argument, we may enunciate 
three plausible conditions for the qq operator in the "maximally attractive 
channel" : 

• qq is gauge invariant 

• qq is a spacetime scalar 

• qq is as local as possible in the lattice x fields 

The <7<72 operator discussed in the previous subsection satisfies each of these 
conditions. For adjoint quarks, however, the Exclusion Principle dictates that 
one of the conditions must be violated, since the equivalent qqa = x tr ( x )x( x ) 
vanishes identically. A possible non-local chirally symmetric condensate could 
form from an operator 

lis =I>„(z)(-l) X * [x tr {x)U»{x) X {x + {i)-x{x)U^x)t r {x + ii)\ . (14) 

This breaks U(2) — > U(l) <8> U(1)a, resulting in two + Goldstones, one of 
which persists even once m, \i 7^ 0. The original U(l)y of (j|) is broken. A 
possible phase diagram resulting from (qq' 3 ) condensation is shown in Fig. 9. 
One interesting issue is that although qq' 3 is a scalar under a lattice version 
of parity,E-3 once a transformation is made to fermion fields with continuum 
degrees of freedom, it is a spacetime pseudoscalar. 

A more compelling possibility is condensation of a gauge non-invariant 
operator 

<?<?3 = X tr ( x ) a ataX( x ), (15) 
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where the t a are generators of the 3, and are antisymmetric in the representa- 
tion in which the are real. Because qq'^ is not gauge singlet, this results in 
a breaking of the SU(2) color group to U(l) by the Higgs mechanism; in other 
words, this is a superconducting solution. A possible phase diagram is shown 
in Fig. 10. Note that since qq^ acts as an adjoint HiggsJield, there is still a 
separation between normal and superconducting phases ,Ej the latter charac- 
terised by a massless photon. By contrast, for QCD in the color-flavor locked 
state, continuity between high and low density phases has been postulatedJla 
We are currently beginning to study the adjoint model using both HMC 
and multi-bosonic methods B3 There are two aspects of the simulation to give 
concern. Firstly, the HMC method permits a minimum Nf = 4, which de- 
stroys asymptotic freedom for adjoint quarks, although that may have little 
direct bearing on the issue of diquark condensation. Secondly, both algo- 
rithms use a power of detM tr M which is guaranteed positive, whereas the 
true measure detM is real but not positive definite. In effect, therefore, the 
simulations incorporate an extra flavor, which allows the possibility of a fla- 
vor non-singlet supcrfluid scalar condensate qq'^ — Xi~ £ ijXj- This condensate 
breaks U(4) — > Sp(4) yielding six Goldstones, some of which must be diquark 
states. Therefore we expect an early onset, and a phase diagram resembling 
Fig. 8. These considerations point to the intriguing possibility of a link be- 
tween the sign problem generic to simulations at fj, ^ and a superconducting 
ground state. 
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